ENUMERATING TOPOLOGICAL (?ife)-CONFIGURATIONS 



JURGEN BOKOWSKI AND VINCENT PILAUD 

Abstract. An (nt;)-configuration is a set of n points and n lines in the projec- 
tive plane such that their point -line incidence graph is fc-regular. The configu- 
ration is geometric, topological, or combinatorial depending on whether lines 
are considered to be straight lines, pseudolines, or just combinatorial lines. 

We provide an algorithm for generating, for given n and k, all topological 
(nj.)-configurations up to combinatorial isomorphism, without enumerating 
first all combinatorial {?ife)-configurations. We apply this algorithm to confirm 
efficiently a former result on topological (184)-configurations, from which we 
obtain a new geometric (184)-configuration. Preliminary results on (194)-confi- 
gurations are also briefly reported. 



1. Introduction 

A point - line configuration is a set P of points and a set L of lines together 
with an incidence relation, where two points of P can be incident with at most 
one line of L and two Unas of L can be incident with at most one point of P. 
Throughout the paper, we only consider connected configurations, where any two 
elements of P U L are connected via a path of incident elements. An isomorphism 
(resp. a duality) between two configurations {P,L) and {P',L') is an incidence- 
preserving map from P Li L to P' U V which sends points to points and lines to 
lines (resp. which exchanges points and lines). 

According to the underlying structure, we distinguish three different levels of 
configurations, in increasing generality: 

Geometric configuration: Points and lines are points and lines in the real pro- 
jective plane P. 

Topological configuration: Points are points in P, but lines are pseudolines, 

i. e. non-separating simple closed curves of P. 
Combinatorial configuration: Just an abstract incidence structure {P, L) as 
described above, with no additional geometric structure. 
In this paper, we focus on regular configurations, i.e. whose incidence relation 
is regular. More precisely, an [uk)- configuration {P,L) is a set P of n points and 
a set L of n lines such that each point of P is contained in k lines of L and each 
line of L contains k points of P. We have represented three famous 3-regular 
configurations in Figure 1 to illustrate the previous definitions. 

Point - line configurations have a long history in discrete 2-dimensional geometry. 
We refer to Branko Griinbaum's recent monograph [Grii09] for a detailed treatment 
of the topic and for historical references. As underlined in this monograph, the 
current study of regular configurations focusses on the following two problems: 
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Figure 1. (Left) Fano's configuration is a combinatorial 
(73)-configuration but is not realizable topologically or geomet- 
rically. (Center) Kantor's topological (103)-configuration is not 
realizable geometrically. (Right) Pappus' configuration is a geo- 
metric (93)-configuration. 



(i) For a given k, determine for which values of n do geometric, topological, and 
combinatorial (nfc)-configurations exist. 

(ii) Enumerate and classify (nfe)-configurations for given k and n. 

In particular, it is challenging to determine the minimal value n for which (nfe)-confi- 
gurations exist and to enumerate these minimal configurations. 

For k S {3,4}, the existence of (nfc)-configurations is almost completely under- 
stood. When k = 3, combinatorial (ri3)-configurations exist for every n> 7, but 
topological and geometric (n3)-configurations exist only for every n > 9. When 
fc = 4, combinatorial (n4)-configurations exist iff n > 13, topological (n4)-configu- 
rations exist iff n > 17 [BS05, BGS09] and geometric (n4)-configurations exist iff 
n > 18 [GriiOe, BSll], with the possible exceptions of 19, 22, 23, 26, 37 and 43. 
For k > 5, the situation is more involved, and the existence of combinatorial, topo- 
logical and geometric (nfc)-configurations is not determined in general. 

Concerning the enumeration, an important effort has been done on combinatorial 
(n^)- and (n4)-configurations. Table 1 provides the known values of the number 
Ck(n) of combinatorial (nj,)-configurations up to isomorphism. The first row of this 
table (k = 3) appeared in [BBPOO], except C3(19) which was announced later on 
in [PBM+04, p.275]. The second row {k = 4) appeared in [BB99, p.34], except 
€4(19) which was only computed recently in [0C12]. 

In this paper, we are interested in the numbers ifc(n) and gkin) of topological 
and geometric (rifc )-configurations up to isomorphism. To obtain these numbers, 
one method is to select the topologically or geometrically realizable configurations 
among the list of all combinatorial (ni:)-configurations. For example, the numbers 
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Table 1. The number Ck{n) of combinatorial (nfc)-configurations 
up to isomorphism. 
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Table 2. The numbers tk{n) of topological (nfc)-configurations 
and gk{n) of geometric (nfe)-configurations up to isomorphism. 



^3(71) and 53 (ji) presented in Table 2 were derived from a careful study of the cor- 
responding combinatorial configurations (see the historical remarks and references 
in [Grii09]). In [Sch07], Lars Schewe provided a general method to study the topo- 
logical realizability of a combinatorial configuration using satisfiability solvers, and 
obtained the numbers t4(17) = 1 and ^4(18) = 16. In [BSll], Jiirgen Bokowski and 
Lars Schewe studied the geometric realizability of a combinatorial configuration. 
This question is clearly an instance of the existential theory of the reals (ETR): it 
boils down to determining whether a set of polynomial equalities and inequalities 
admits a solution in the reals (indeed, the inclusion of a point in a line can be tested 
by a polynomial equation). Using the construction sequences presented in [BSll], 
the complexity of this instance of ETR can be decreased significantly. With this 
method, Jiirgen Bokowski and Lars Schewe showed that the only combinatorial 
(174)-configuration which is topologically realizable is not geometrically realizable 
and they exhibited a geometric (184)-configuration. 

Table 2 summarizes the values of t^{n), g^(ri),ti(ri) and gi{n) known up-to-date 
(we have additionally included our results in bold letters; see below). This table 
indicates a clear difference of behavior between 3- and 4-regular configurations. 
On the one hand, when A; = 3, most of the combinatorial (n3)-configurations are 
topologically and geometrically realizable for small values of n. For n < 12, the only 
counter-examples are the Fano (73)-configuration, the Mobius-Kantor (83)-configu- 
ration, and Kantor's (103)-configuration — see Figure 1 (left & center). On the 
other hand, when fc = 4, it is not reasonable to look for geometric (n4)-configura- 
tions among all combinatorial (n4)-configurations. To further extend our knowledge 
on geometric configurations, it thus seems crucial to limit our research to those 
combinatorial configurations which are already topologically realizable. 

Motivated by this observation, we present an algorithm for generating, for given n 
and fc, all topological (nfc)-configurations up to isomorphism, without enumerating 
first all combinatorial (nfe)-configurations. The algorithm sweeps the projective 
plane to construct a topological (nfc)-configuration {P,L), but only considers as 
relevant the events corresponding to the sweep of points of P. This strategy enables 
us to identify along the way some isomorphic topological configurations, and thus 
to maintain a reasonable computation space and time. 

We have developed two different implementations of this algorithm. The first 
one was written in HASKELL by the first author to develop the strategy of the 
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enumeration process. Once the general idea of the algorithm was settled, the second 
author wrote another implementation in JAVA, focusing on the optimization of 
computation space and time of the process. 

We outline three applications of our algorithm. First, the algorithm is interesting 
in its own right. Before describing some special methods for constructing topological 
configurations, Branko Griinbaum writes in [Gru09, p. 165] that "the examples of 
topological configurations presented so far have been ad hoc, obtained essentially 
through (lots of) trial and error". Our algorithm can reduce considerably the "trial 
and error" method. Second, our algorithm enables us to check and confirm all 
values of i4(n), for n < 18, obtained in earlier papers. We can use for that a single 
method and reduce considerably the computation time {e.g. the computation of the 
(184)-configurations needed several months of CPU-time in [Sch07], and only one 
hour with our JAVA implementation). Finally, this algorithm enables us to compute 
all ^4(19) — 4028 isomorphism classes of topological (194)-configurations. 

As an application of our enumeration results, we studied in detail the possible 
geometric realizations of the topological (184)-configurations. Using a maple in- 
plementation of the construction sequence method of Jiirgen Bokowski and Lars 
Schewe [BS05], we obtain that there are precisely 2 geometric (184)-configurations: 
the first (184)-configuration constructed in [BS05], plus an additional one which 
appears for the first time in this paper. In contrast, deriving the list of geometric 
(194)-configurations from the list of topological (194)-configurations still requires 
some computational effort and is left to a subsequent paper. 

The first section of this paper is devoted to the enumeration algorithm for iso- 
morphism classes of topological configurations. The second section presents the 
application to the enumeration of geometric (184)-configurations. 

Topological configurations are pseudoline arrangements, or rank 3 oriented ma- 
troids. We assume the reader to have some basic knowledge on these topics. We 
refer to [Bok06, BLS+99, Knu92] for introductions. 

2. Topological configurations 

In this section, we present our algorithm to generate all isomorphism classes of 
topological (nfc)-configurations, for given n and k. Let us insist again on the crucial 
fact that we do not need to enumerate first all combinatorial (nfe)-configurations. 
The main idea of the algorithm is to sweep the projective plane to construct a 
topological (nfc)-configuration {P,L), only focussing on the "relative positions of 
the points of P" and ignoring at first the "relative positions of the other crossings 
of the pseudolines of L" (precise definitions are given in Section 2.1). This strategy 
enables us to identify along the way some isomorphic topological configurations, 
and thus to maintain a reasonable computation space and time. 

2.1. Three equivalence relations. There are three distinct notions of equiva- 
lence on topological configurations. 

The finest notion is the usual notion of topological equivalence between pseu- 
doline arrangements in the projective plane: two configurations are topologically 
equivalent if there is an homeomorphism of their underlying projective planes that 
sends one arrangement onto the other. 

The coarsest notion is that of combinatorial equivalence: two (nfc)-configurations 
are combinatorially equivalent if they are isomorphic as combinatorial (nfc)-confi- 
gurations. 
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The intermediate notion is based on the graph of admissible mutations. Re- 
member that a mutation in a pseudohne arrangement is a local transformation of 
the arrangement where only one pseudoline i moves, sweeping a single vertex v of 
the remaining arrangement. It only changes the position of the crossings of i with 
the pseudolines incident to v. If those crossings are all 2-crossings, the mutation 
does not perturb the /c-crossings of the arrangement, and thus produces another 
topological (nfc)-configuration. We say that such a mutation is admissible. Two 
configurations are mutation equivalent if one can be obtained from the other by a 
(possibly empty) sequence of admissible mutations followed by an homcomorphism 
of the underlying projective space. 




Figure 2. An admissible mutation. 

Obviously, topological equivalence implies mutation equivalence, which in turn 
implies combinatorial equivalence. The reciprocal implications are wrong. 

As an illustration, the two (184)-configurations depicted in Figure 3 are com- 
binatorially equivalent (the labels on the pseudolines provide a combinatorial iso- 
morphism) but not topologically equivalent (the left one has 22 quadrangles and 
2 pentagons, while the right one has 23 quadrangles). In fact, one can even check 
that they are not mutation equivalent. 




Figure 3. Two (184)-configurations which are combinatorially 
equivalent but neither mutation nor topologically equivalent. 
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2.2. Representation of arrangements. In this section, we state certain proper- 
ties of configurations that we can assume without loss of generahty. In particular, 
we choose a suitable representation of our pseudoline arrangements that we will use 
for the description of the algorithm in Section 2.3. 

Simple configurations — A topological configuration (P, L) is simple if no 
three pseudolines of L meet at a common point except if it is a point of P. Since 
any topological (n/j)-configuration can be arbitrarily perturbed to become simple, 
we only consider simple topological (rifc )-configurations. Once we obtain all simple 
topological (n^.)-configurations, it is usual to obtain all (non- necessarily simple) 
topological (n;j)-configurations up to topological equivalence by exploring the mu- 
tation graph, and we do not report on this aspect. 

In a simple (nfe)-configuration (P, i), there are two kinds of intersection points 
among pseudolines of L: the points of P, which we also call k-crossings, and the 
other intersection points, which we call 2-crossings. Each pseudoline of L contains 
k fc-crossings and n — 1 — k{k — 1) 2-crossings. In total, a simple (7i;j)-configuration 
has n fc-crossings and (j) — ^((2) ^ 1) 2-crossings. 

Segment length distributions — A segment of a topological configuration (P, L) 
is the portion of a pseudoline of L located between two consecutive points of P. 
If (P, L) is simple, a segment contains no fc-crossing except its endpoints, but may 
contain some 2-crossings. The length of a segment is the number of 2-crossings it 
contains. 

The circular sequence of the segment lengths on a pseudoline of L forms a 
fc-partition of rt — 1 — fc(fc — 1). We call a maximal representative of a fc-tuplc the 
lexicographic maximum of its orbit under the action of the dihedral group {i.e. rota- 
tions and refiections of the fc-tuple) . We denote by 11 the list of all distinct maximal 
representatives of the fc-partitions of n — 1 — k{k — 1), ordered lexicograhically. For 
example, when fc = 4 and n = 17, we have H = [4,0,0,0], [3,1,0,0], [3,0,1,0], 
[2,2,0,0], [2,0,2,0], [2,1,1,0], [2,1,0,1], [1,1,1,1]. 

A SUITABLE representation — We represent the projective plane as a disk where 
we identify antipodal boundary points. Given a simple topological (rifc)-configura- 
tion (P, L), we fix a representation of its underlying projective plane which satisfies 
the following properties (see Figure 4 left). 

The leftmost point of the disk (which is identified with the rightmost point of the 
disk) is a point of P, which we call the base point. The k pseudolines of L passing 
through the base point are called the frame pseudolines., while the other n — k 
pseudolines of L are called working pseudolines. The frame pseudolines decompose 
the projective plane into fc connected regions which we call frame regions. A crossing 
is a frame crossing if it involves a frame pseudoline and a working crossing if it 
involves only working pseudolines. 

The boundary of the disk is a frame pseudoline, which we call the base line. 
We furthermore assume that the segment length distribution A on the top half- 
circle appears in 11 {i.e. is its own maximal representative), and that no maximal 
representative of the segment length distribution of a pseudoline of L appears before 
A in n. In particular, the leftmost segment of the base line is a longest segment of 
the configuration. 
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Figure 4. Suitable representation of a (174)-configuration, and 
the corresponding wiring diagram. 
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Wiring diagram and allowable sequence — Another interesting representa- 
tion of our (rifc )-configuration is the wiring diagram [GP93] of its working pseudo- 
lines (see Figure 4 right) . It is obtained by sending the base point to infinity in the 
horizontal direction. The frame pseudolines are k horizontal lines, and the n — k 
working pseudolines are vertical wires. The orders of the working pseudolines on a 
horizontal line sweeping the wiring diagram from top to bottom form the so-called 
allowable sequence of the working arrangement, as defined in [GP93]. 

2.3. Description of the algorithm. Our algorithm can enumerate all topological 
(nfc)-configurations up to either topological or combinatorial equivalence. In order 
to maintain a reasonable computation space and time, the main idea is to focus on 
the relative positions of the points of the configurations and to ignore at first the 
relative positions of the other crossings among the pseudolines. In other words, to 
work modulo mutation equivalence as defined in Section 2.1. 

More precisely, we first enumerate at least one representative of each mutation 
equivalence class of topological (nfe)-configurations. From these representatives, we 
can derive: 

(1) all topological (rifc)-configurations up to topological equivalence: we explore 
each connected component of the mutation graph with our representatives 
as starting nodes. 

(2) all combinatorial (7ife)-configurations that are topologically realizable: we 
reduce the result modulo combinatorial equivalence. 

Since our motivation is to study geometric (nfc)-configurations, we are only inter- 
ested by point (2). We discuss a relatively efficient approach to test combinatorial 
equivalence in Section 2.4. In this section, we give details on the different steps in 
our algorithm. 
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Sweeping process — Our algorithm sweeps the projective plane to construct a 
topological (nfc)-configuration. The sweep line sweeps the configuration from the 
base line on the top of the disk to the base line on the bottom of the disk. Inside 
each frame region, it always passes through the base point and always completes 
the configuration into an arrangement of n + 1 pseudolines. When it switches from 
one frame region to the next one, it coincides with the separating frame pseudoline. 
Along the way, it sweeps completely all the working pseudolines. Except those 
located on the frame pseudolines, we assume that the crossings of the configuration 
are reached one after the other by the sweep line. After the sweep line swept a 
crossing, we remember the order of its intersections with the working pseudolines. 
In other words, the sweeping process provides us with the allowable sequence of the 
working pseudolines of our configuration. 

Since admissible mutations are irrelevant for us, we only focus on the moments 
when our sweep line sweeps a fc-crossing. Thus, two different events can occur: 

• when the sweep line sweeps a working A;-crossing, and 

• when the sweep line sweeps a frame pseudoline. 

In the later case, we sweep simultaneously k — 1 frame fc-crossings (each involving 
the frame pseudoline and k — 1 working pseudolines), and n — 1 — fc(fc — 1) frame 
2-crossings (each involving the frame pseudoline and a working pseudoline). Be- 
tween two such events, the sweep line may sweep working 2-crossings which are 
only taken into account when we reach a new event. Let us repeat again that the 
precise positions of these working 2-crossings is irrelevant in our enumeration. 

To obtain all possible solutions, we maintain a stack with all subconfigurations 
which have been constructed so far, remembering for each one: 

(i) the order of the working pseudolines on the current sweep line, 

(ii) the number of frame and working fc-crossings and 2-crossings which have al- 
ready been swept on each working pseudoline, 

(iii) the length of the segment currently swept by the sweep line, and 

(iv) the history of the sweeps performed to reach this subconfiguration. 

At each step, we remove the first subconfiguration from the stack, and insert all ad- 
missible subconfigurations which can arise after sweeping a new working fc-crossing 
or a new frame pseudoline. We finally accept a configuration once we have swept 
fc frame pseudolines and n — k{k — 1) — 1 working fc-crossings. 

Any subconfiguration considered during the algorithm is a potential (nfe)-confi- 
guration. Throughout the process, we make sure that any pair of working pseudo- 
lines cross at most once, that the number of frame pseudolines (resp. of working 
fc-crossings) already swept never exceeds fc (resp. n — 1 — k{k — 1)), and that the 
total number of working 2-crossings never exceeds in — 2k){n — 1 — fc(fc — l))/2. 
Furthermore, on each pseudoline, the number of frame and working fc-crossings 
(resp. 2-crossings) already swept never exceeds fc (resp. n — 1 — fc(fc — 1)), the num- 
ber of working 2- and fc-crossings already swept never exceeds n — I — k{k — 1), and 
the segment currently swept is not longer than the leftmost segment of the base line. 

We now detail individually each step of the algorithm. 

Initialization — We initialize our algorithm sweeping the base line. We only 
have to choose the distribution of the lengths of the segments on the base line. The 
possibilities are given by the list 11 of maximal representatives of fc-partitions of 
n - 1 - fc(fc - 1). 
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Figure 5. Sweeping a working fc-crossing. 

Sweep a working /c-crossing — If we decide to sweep a working fc-crossing, 
we have to choose the k working pseudolines which intersect at this fc-crossing, and 
the direction of the other working pseudolines. 

Since we are allowed to perform any admissible mutation, we can assume that 
all the pseudolines located to the left of the leftmost pseudoline of the working 
fc-crossing, and all those located to the right of the rightmost pseudoline of the 
working fc-crossing do not move. 

We say that the pseudolines located between the leftmost and the rightmost 
pseudolines of the working fc-crossing form the kernel of the working fc-crossing. 
We have to choose the positions of the pseudolines of the kernel after the flip: each 
pseudoline of the kernel either belongs to the working fc-crossing, or goes to its left, 
or goes to its right (see Figure 5). 

A choice of directions for the kernel is admissible provided that 

(i) each pseudoline involved in the fc-crossing can still accept a working fc-crossing; 

(ii) each pseudoline of the kernel can still accept as many working 2-crossings as 
implied by the choice of directions for the kernel; 

(iii) no segment becomes longer than the leftmost segment of the base line; and 

(iv) any two pseudolines which are forced to cross by the choice of directions for 
the kernel did not cross earlier {i.e. they still form an inversion on the sweep 
line before we sweep the working fc-crossing). 

Sweep a frame pseudoline — If we decide to sweep a frame pseudoline, we 
have to choose the (fc — 1)^ working pseudolines involved in one of the fc — 1 frame 
fc-crossings, and the direction of the other working pseudolines. 

As before, we can assume that a pseudoline does not move if it is located to the 
left of the leftmost pseudoline involved in one of the fc — 1 frame fc-crossings, or to 
the right of the rightmost pseudoline involved in one of the fc — 1 frame fc-crossings. 
Otherwise, we can perform admissible mutations to ensure this situation. 

The other pseudolines form again the kernel of the frame sweep, and we have to 
choose their positions after the flip. Each pseudoline of the kernel either belongs 
to one of the fc — 1 frame fc-crossings, or can choose among fc possible directions: 
before the first frame fc-crossing, or between two consecutive frame fc-crossings, or 
after the last frame fc-crossing (see Figure 6). 

As before, a choice of directions for the kernel is admissible if 

(i) each pseudoline involved (resp. not involved) in one of the fc — 1 frame 
fc-crossings can still accept a frame fc-crossing (resp. a frame 2-crossing); 

(ii) each pseudoline of the kernel can still accept as many working 2-crossings as 
implied by the choice of directions for the kernel; 
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Figure 6. Sweeping a frame pseudoline (right). 

(iii) no segment becomes longer than the leftmost segment of the base line; and 

(iv) any two pseudolines which are forced to cross by the choice of directions for 
the kernel did not cross earlier {i.e. they still form an inversion on the sweep 
line before we sweep the frame pseudoline). 

Sweep the last frame region — Our sweeping process finishes once we have 
swept n — 1 — k(k — 1) working fc-crossings and k frame pseudolines. Each resulting 
subconfiguration should still be completed into a topological (n/j)-configuration 
with some necessary remaining 2-crossings. More precisely, we need to add on 
each working pseudoline as many working 2-crossings as its number of inversions in 
the permutation given by the working pseudolines on the final sweep line, without 
creating segments that are too long. 

After this last selection, all the constructed configurations are finally guaranteed 
to be valid topological (nfc)-configurations. To make sure that we indeed obtain 
the representation presented in Section 2.2, we remove each configuration (P, L) in 
which the maximal representative of the segment length distribution of a pseudoline 
of L appears in the list 11 before the segment length distribution of its base line. 

2.4. Testing combinatorial equivalence. In Section 2.1, we have seen three 
equivalence relations between topological (nfc)-configurations: combinatorial, muta- 
tion and topological equivalence. As explained in Section 2.3, our algorithm outputs 
at least one representative per mutation equivalence class of topological (nfe)-confi- 
gurations. However, we can obtain more than one representative per class, and 
two topological (rifc)-configurations which are not mutation equivalent can still be 
combinatorially equivalent. We thus need to reduce the output of our algorithm. 

Note that the topological equivalence between two (nfe)-configurations (Pi,Li) 
and {P2,L2) can be tested in Q{n^) time. Indeed, since the topological configura- 
tions are embedded on the projective plane, the matchings between Pi and P2 and 
between Li and L2 induced by an homeomorphism mapping (Pi, Li) to {P2, L2) are 
determined by the images of any two distinguished pseudolines of Li. There- 
fore, for each of the 8(71^) possible choices for the images of i^l' , we can test in 
linear time whether this choice yields or not an homeomorphism between (Pi,Li) 
and (P2, ^2)- Both combinatorial and mutation equivalences are however harder to 
decide computationally. We focus here on methods and heuristics to quickly test 
combinatorial equivalence. 

In order to limit unnecessary computation, we make use of combinatorial invari- 
ants associated to configurations. If two configurations have distinct invariants, 
they cannot be combinatorially equivalent. Reciprocally, if they share the same in- 
variant, it provides us with information on the possible combinatorial isomorphisms 
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between these two configurations. The invariants we have chosen are the clique and 
coclique distributions. We furthermore need a multiscale invariant technique, based 
on the notion of derivation of a combinatorial invariant. We introduce these notions 
and methods in the next paragraphs. 

Cliques and COCLIQUES — Let (P, £) be a combinatorial configuration. For 
j > 3, define a j-clique of (P, L) to be any set of j points of P which are pairwise 
related by lines of L. For any point p of P, let 7j(p) be the number of j-cliques 
containing p, and let 7(p) := {'lj{p))j>3- The clique distribution of {P, L) is the 
multiset -f{P) lj{p) \ p G P}. 

Similarly, a j -coclique of {P, L) is a set of j lines of L which are pairwise inter- 
secting at points of P. For any line £ of L, let 6j{£) be the number of j-cocliques 
containing £, and let 6{£) :— {Sj{£))j>3. The coclique distribution of {P, L) is the 
multiset d{L) :— ^S{£) \ £ e L^. In other words, the coclique distribution of {P,L) 
is the clique distribution of its dual configuration (L,P). 

The pair (7(P), S{L)) of clique and coclique distributions of the configuration {P, L) 
is a natural and powerful combinatorial invariant of (P, L). 

Derivation of combinatorial invariants — Let (P, L) be a (7T,fc)-configura- 
tion. Assume that 7 : P — >■ X and 5 : L ^ Y are two functions from the point 
set and the line set of (P, L) respectively to arbitrary sets X and F, such that 
the multisets 7(P) := {{7(p) | p G P^- C X and 5{L) := 15{£) \ £ d Y sue 

combinatorial invariants of (P, L). The clique and coclique distributions are typical 
examples of such functions 7 and 5. Observe that we again abuse notation: the 
functions 7 and 5 usually depend on the configuration (P, L), but we consider that 
this dependence is clear from the context. Note however that the target sets X 
and y of 7 and 5 do not depend upon (P, L). 

While reducing a set of configurations up to combinatorial equivalence, such a 
pair of combinatorial invariants (7(P), 5{L)) can be used in two different ways: 

(i) either to separate classes of combinatorial isomorphism: two configurations 
with different invariants cannot be combinatorially equivalent; 

(ii) or to guess combinatorial isomorphisms: an isomorphism between two config- 
urations should respect the invariants 7 and 5. 

It often happens however that the pair of combinatorial invariants {'y{P),S{L)) is 
not precise enough neither to distinguish two configurations, nor to guess a com- 
binatorial isomorphism between them. It occurs when many points (resp. many 
lines) of a configuration (P, L) get the same image under 7 (resp. under i5). Two 
fundamentally different cases can lead to this situation. On the one hand, the con- 
figuration (P, L) can have a large automorphism group. In this case, points (resp. 
lines) in a common orbit under the automorphism group cannot be distinguished 
combinatorially, and thus no invariant can speed up the isomorphism test. On the 
other hand, it could also be that the combinatorial invariant {'j{P),5{L)) is not 
precise enough to distinguish the neighborhood properties of the points with the 
same image under 7 (resp. the lines with the same image under S). In the later case, 
we can construct a new pair of combinatorial invariants which refines (7(P), S{L)), 
taking into account the neighborhoods of points and lines in the configuration. We 
call these invariants the derivatives of 7 and S and denote them 7' and S' . 

The derivative of the invariant 7 : P — >■ X is the function 7' : L — >■ X'^ which 
associates to a line t oi L the multiset ^'{£) '■— \l{p) \ p ^ P,p ^ £^- Intuitively, 
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the image "/'{£) of a line £ contains all the combinatorial information carried by 7 
concerning the points of P contained in £. Similarly, the derivative of the invariant 
5 : L —i' Y is the function 5' : P which associates to a point poiP the multiset 

5'{p) := 15{£:) \£eL,pe £}. The pair {S' (P) , (L)) is a pair of combinatorial 
invariants as defined previously, and it refines the previous pair {j{P), 6{L)). 

If this new invariant is still not precise enough, we can consider higher order 
derivatives 7^") := (7("-i))' and J^") := (J^""!))' of the initial invariants. We 
obtain this way a family of refinements of {j{P), S{L)). Of course, these invariants 
ultimately carry the same combinatorial information. We use this family in the 
following multiscale technique. 

MuLTiSCALE INVARIANTS — The main idea of our reduction process is to use 
derivative invariants in a multiscale process. Consider a set C of configurations that 
we want to reduce up to combinatorial equivalence. Assume that 7 : P — > X and 
S : L ^ Y are two functions defining a pair of combinatorial invariants {"/{P), S{L)) 
of a configuration (P, L) . We separate the configurations of C into classes with 
distinct invariants, which we can consider independently. We now compute the 
derivative invariants ((5'(P), 7'(L)) for each configuration {P,L). For a given class, 
we then have three possible situations: 

(1) If the derivative invariants {6' (P) , -f' (L)) are not the same for all config- 
urations (P, L) of the class, we split the class into refined subclasses and 
reiterate the refinement (computing one more derivative). 

(2) If the derivative invariants ((5'(P), 7'(L)) are the same for all configura- 
tions (P, L) of the class but determine more information on the possible 
isomorphisms between configurations of the class than the original invari- 
ants (7(P),5(L)), then we reiterate the refinement. 

(3) Otherwise, the derivative invariants (5'(P), 7'(L)), as well as any further 
derivative, provide the same combinatorial information as the original in- 
variants (7(P), (5(L)). Thus, we stop the refinement process and start a 
brute-force search for possible isomorphisms between the remaining con- 
figurations in the class. The efficiency of this brute-force search depends 
on the quality of the combinatorial information provided by the invari- 
ants (7(P),5(i)). 

This process can be seen as a multiscale process: typically, some invariants provide 
sufficiently information to deal with certain classes of C, while other classes require 
far more precision (obtained by derivatives) to be reduced. 

Using this multiscale technique, starting from the clique and cocliquc distribu- 
tions of configurations, we managed to reduce the 69 991 topological (194)-confi- 
gurations produced by our sweeping algorithm into 4 028 classes of combinatorial 
equivalence in about one hour^. 

2.5. Results. We present in this section the results of our algorithm. First, it 
enables us to check efficiently all former enumerations of topological (nfc)-confi- 
gurations. The JAVA implementation developed by the second author finds all 
(nfc)-configurations in less than a minute^ when k = 3 and n < 11, or when k — 4 
and n < 17. In particular, we checked that there is no topological (7i4)-configura- 
tion when n < 16 [BS05], and that there is a single topological (174)-configuration 



^Computation times on a 2.4 GHz Intel Core 2 Duo processor with 4Go of RAM. 
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up to combinatorial isomorphism. This configuration is represented in Figure 7, 
and labeled in such a way that: 

• the quarter-turn rotation which generates the symmetry group of the pic- 
ture is the permutation (A)(B,C,D,E)(F,G,H,I)(J,K,L,M)(N,0)(P,Q); and 

• the permutation (A,a)(B,b) . . . (P,p)(Q,q) is a self-polarity of the topolog- 
ical configuration. 
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Figure 7. The topological (174)-configuration [BGS09]. 



When k = A and n — 18, we reconstructed the 16 combinatorial equivalence 
classes of topological (184)-configurations obtained in [Sch07] with satisfiability 
solvers. See [BSll, Figure 6] for a description of these configurations. To obtain this 
result, our implementation needed about one hour^, compared to several months 
of CPU-time required in [Sch07]. The two (184)-configurations presented in Fig- 
ure 3, which are combinatorially equivalent but not mutation equivalent, occurred 
while we were reducing the list of (184)-configurations up to combinatorial equiva- 
lence, using as a first reduction a certain invariant of mutation equivalence defined 
in [BS12]. In the next section, we present two combinatorially distinct geometric 
(184)-configurations obtained from the list of topological (184)-configurations. 

Finally, we want to report on preliminary results concerning the enumeration of 
topological (194)-configurations, which initially motivated our work. In about 15 
days of computation time^, we obtained the complete list of topological (194)-confi- 
gurations: 



^Computation times on a 2.4 GHz Intel Core 2 Duo processor with 4Go of RAM. 
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Result 1. There are precisely 4028 topological {19 4) -configurations up to combi- 
natorial equivalence. Among them, 222 are self-dual. 

From this list, we can immediately extract examples of topological (1 94) -configu- 
rations with non-trivial symmetry groups, closing along the way an open question 
of Branko Griinbaum [Gru09, p. 169, Question 5]. The next step is naturally to 
study the possible geometric realizations of all these topological (194)-configura- 
tions. This work in progress still requires an important computational effort and 
will be reported in a subsequent paper. 

3. Application to geometric (184)-configurations 

As an application of the enumeration of topological configurations, we derive 
all isomorphism classes of geometric (184)-configurations. To obtain it, we im- 
plemented in MAPLE the construction sequence method of Jiirgen Bokowski and 
Lars Schewe [BSll]. Among the 16 topological (184)-configurations (first gener- 
ated by Lars Schewe [Sch07] and now confirmed by our JAVA program), only 8 are 
compatible with Pappus' and Desargues' Theorem. Starting from these remaining 
configurations, we run our MAPLE code and obtain the following result: 

Result 2. There are precisely two geometric (IS 4) -configurations up to combina- 
torial isomorphism. 

The first geometric (184)-configuration was obtained in [BS05, Section 4]. 
In Figure 8, we have labeled its points A, . . . , R and lines a, . . . , r in such a way 
that the permutation (A,a)(B,b) . . . (Q,q)(R,r) is a self-duality of the configuration. 
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Figure 8. Bokowski and Schewe's geometric (184)-configuration [BS05]. 
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The automorphism group of the combinatorial configuration is generated by the 
permutations: 

(A,B,C)(D,E,F)(G,H,I)(J,K,L)(M,N,0)(P,Q,R) 
(A)(K)(B,C,L,J)(D,F,I,R)(E,Q,M,G)(H,0,N,P) 
(A)(K)(B,L)(C,J)(D,I)(E,M)(F,R)(G,Q)(H,N)(0,P) 
and is isomorphic to the symmetric group on 4 elements. Together with the self- 
duality (A,a)(B,b) . . . (Q,q)(R,r), the automorphism group of the Levi graph of the 
configuration is thus isomorphic to 64 x Z2. Observe that only the first permuta- 
tion (A,B,C) . . . (P,Q,R) and the duality (A,a)(B,b) . . . (Q,q)(R,r) are geometrically 
visible, while the other generators of the automorphism group of the combinato- 
rial configuration are not isometrics of the geometric configuration of Figure 8. In 
Figure 9, we have performed a projective transformation of the configuration of 
Figure 8 (sending the four 3-valent points in Figure 8 to a square). The last gener- 
ator (A)(K)(B,L) . . . (0,P) then becomes a central symmetry in the new geometric 
(184)-configuration. 

The realization space of this configuration consists of two points, both expressed 
with coordinates in Q [l + a/5] . 





Figure 9. Another geometric realization of Bokowski and 
Schewe's geometric (184)-configuration [BS05] of Figure 8. 

The second geometric (184)-configuration is a result of our maple code and 
appears for the first time in this paper. 

In Figure 10, we have labeled its points A, . . . , R and lines a, . . . , r in such a way 
that the permutation (A,a)(B,b) . . . (Q,q)(R,r) is a self-polarity of the configuration. 

The automorphism group of the combinatorial configuration is generated by the 
permutation (Q)(R)(A,P)(B,0)(C,N)(D,M)(E,L)(F,K)(G,J)(H,I). Together with the 
self-polarity (A,a)(B,b) . . . (Q,q)(R,r), the automorphism group of the Levi graph 
of the configuration is thus isomorphic to Z2 x Z2 . This group is completely realized 
in the geometric representation of Figure 10. 

The realization space of this configuration consists of two points, both expressed 



with coordinates in Q v^l08 -I- 12^93 
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Figure 10. The new geometric (184)-configuration. 



To conclude, we want to emphasize that the discovery of the first (184)-configu- 
ration of Figure 8 inspired Branko Griinbaum to find a new family of (6rn4)-confi- 
gurations, for any m > 3 (see [Grii09, Chapter 3, p. 171] and Figure 11). This 
raises the following appealing open question: 

Problem 3. Generalize our second geometric {I84) -configuration of Figure 10 to 
obtain another new infinite family of geometric [n 4) -configurations. 




Figure 11. The (6m)-family inspired by the geometric 
(184)-configuration of Figure 8. 
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For example, we have been able to derive from the second geometric (184)-confi- 
guration of Figure 10 a family of ((18 + 17m)4)-configurations. Unfortunately, the 
set 18 + 17N does not intersect the set {19, 22, 23, 26, 37, 43} of values n for which 
no (n4)-configuration is known. 

Acknowledgements 

The first author thanks three colleagues from the Universidad Nacional Autonoma 
de Mexico, namely Ricardo Strausz Santiago, Rodolfo San Augustin Chi, and Oc- 
tavio Paez Osuna, for many stimulating discussions about various different earlier 
versions of the presented algorithm during his one year sabbatical stay (2008/2009) 
in Mexico City. We also thank Leah Herman from the University of Alaska Fair- 
banks for valuable discussions and comments about the subject. We are grateful to 
Branko Griinbaum, Tomaz Pisanski, and Gunnar Brinkmann for encouragements 
and helpful communications. As frequent users, we are indebted to the development 
team of the geometric software Cinderella, in particular Jiirgen Richter-Gebert 
and Ulrich Kortenkamp. Finally, we thank two anonymous referees for their com- 
ments and suggestions on the presentation. 

References 

[BB99] Anton Betten and Dieter Betten. Tactical decompositions and some configurations V4. 
J. Geom., 66(l-2):27-41, 1999. 

[BBPOO] Anton Betten, Gunnar Brinkmann, and Tomaz Pisanski. Counting symmetric con- 
figurations 1)3. Discrete Appl. Math., 99(1-3):331— 338, 2000. Proceedings of the 5th 
Twente Workshop on Graphs and Gombinatorial Optimization (Enschede, 1997). 

[BGS09] Jiirgen Bokowski, Branko Griinbaum, and Lars Scfiewe. Topological configurations 
{n4) exist for all n > 17. European J. Gombin., 30(8): 1778-1785, 2009. 

[BLS+99] Anders Bjorner, Micliel Las Vergnas, Bernd Sturmfels, Neil White, and Giinter M. 

Ziegler. Oriented matroids, volume 46 of Encyclopedia of Mathematics and its Appli- 
cations. Cambridge University Press, Cambridge, second edition, 1999. 

[Bok06] Jiirgen Bokowski. Gomputational oriented matroids. Cambridge University Press, 
Cambridge, 2006. 

[BS05] Jiirgen Bokowski and Lars Schewe. There are no realizable I54- and 164-configurations. 
Rev. Roumaine Math. Pures Appl, 50(5-6):483-493, 2005. 

[BSll] Jiirgen Bokowski and Lars Schewe. On the finite set of missing geometric configurations 
(714). To appear in Gomputational Geometry: Theory and Applications, 2011. 

[BS12] Jiirgen Bokowski and Ricardo Strausz Santiago. A manifold associated to a topological 
(nfe)-configuration. Preprint, 2012. 

[GP93] Jacob E. Goodman and Richard Pollack. Allowable sequences and order types in dis- 
crete and computational geometry. In New trends in discrete and computational ge- 
ometry, volume 10 of Algorithms Gombin., pages 103-134. Springer, Berlin, 1993. 

[Grii06] Branko Griinbaum. Connected (714) configurations exist for almost all n — second up- 
date. Geombinatorics, 16(2):254-261, 2006. 

[Grii09] Branko Griinbaum. Gonfigurations of points and lines, volume 103 of Graduate Studies 
in Mathematics. American Mathematical Society, Providence, RI, 2009. 

[Knu92] Donald E. Knuth. Axioms and hulls, volume 606 of Lecture Notes in Gomputer Sci- 
ence. Springer- Verlag, Berlin, 1992. 

[OC12] Octavio Paez Osuna and Rodolfo San Agustin Chi. The combinatorial (I94) configu- 
rations. Ars Math. Gontemp., 5(2):231-237, 2012. 

[PBM+04] Tomaz Pisanski, Marko Boben, Dragan Marusic, Alen Orbanic, and Ante Graovac. 
The 10-cages and derived configurations. Discrete Math., 275(1-3) :265-276, 2004. 

[Sch07] Lars Schewe. Satisfiability Problems in Discrete Geometry. PhD thesis, Technische 
Universitat Darmstadt, 2007. 



18 



JURGEN BOKOWSKI AND VINCENT PILAUD 



Tbchnische Universitat Darmstadt 

E-mail address: juergen.bokowskiagooglemail.com 

CNRS & LIX, ECOLE POLYTECHNIQDE, PaLAISEAD 
E-mail address: vincent.pilaudOlix.polytechnique.fr 
URL: http : //www . lix . polytechnique . f r/"pilaud/ 



